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Suppressing longitudinal double-layer oscillations by using elliptically polarized 
laser pulses in the hole-boring radiation pressure acceleration regime 
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It is shown that well collimated mono-energetic ion beams with a large particle number can be generated 
in the hole-boring radiation pressure acceleration regime by using an elliptically polarized laser pulse with 
appropriate theoretically determined laser polarization ratio. Due to the J x B effect, the double-layer charge 
separation region is imbued with hot electrons that prevent ion pileup, thus suppressing the double-layer 
oscillations. The proposed mechanism is well confirmed by Particle-in-Cell simulations, and after suppressing 
the longitudinal double-layer oscillations, the ion beams driven by the elliptically polarized lasers own much 
better energy spectrum than those by circularly polarized lasers. 

PACS numbers: 52.38.Kd, 41.75.Jv, 52.35.Mw, 52.59.-f 



Well-collimated mono-energetic ion beams are useful 
for producing high energy density matter, radiograph- 
ing transient processes, tumor therapy, and ion-fast ig- 
nition in laser fusion^—. Radiation pressure acceleration 
(RPA) is an efficient scheme for generating high quality 
ion beams, which as a particular case of laser-induced 
cavity pressure acceleration (LICPA) 7 ~ 9 usually call for 
relatively high laser intensity. According to the target 
thickness, usually there are two modes of RPA accelera- 
tion mechanisms: light sail (LS) RPA for thin target^— 
and hole boring (HB) RPA for thick target^—. In par- 
ticular, the HBRPA owns the intrinsic property for large 
particle number acceleration^!. In HBRPA, the pondero- 
motive force of a laser pulse drives most or all the local 
electrons inward, resulting in a shock-like double-layer 
(DL) region with huge electrostatic charge-separation 
field. The latter traps and reflects the ions ahead of 
the DL, and compresses and accelerates them inward like 
a piston^r— . For RPA usually a circularly polarized 
(CP) laser pulse is invoked^—, since its ponderomotivc 
force does not have a high-frequency oscillating compo- 
nent that can preheat the electrons and thus reduce the 
electrostatic DL field and efficiency of ion acceleration. 
However, as pointed out by previous works 20 i 21 i 25 , the 
intense electrostatic DL field oscillates because trapping 
and reflection of ions by the this field are accompanied 
by consecutive increase and decrease of the ion therein. 
The oscillations lead to significant increase of the energy 
spread of the accelerated ions. On the other hand, if 
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a linearly polarized (LP) laser pulse is used, the oscil- 
lating electrostatic DL field will be smeared out by the 
hot electrons generated by the second-harmonic compo- 
nent of the ponderomotive force. However, the DL struc- 
ture is also destroyed and no mono-energetic ion beam is 
generated 2 ^. 

In this paper, we propose to use an elliptically polar- 
ized (EP) laser pulse for ion acceleration in the HBRPA 
regime. In this scheme, J x B electron heating causes 
the DL region to be filled with an appropriate distri- 
bution of hot electrons such that large ion density vari- 
ations and electrostatic field oscillations in the DL re- 
gion are suppressed, so that a mono-energetic ion beam 
can be generated. An analytical model is used to esti- 
mate the most suitable laser polarization ratio needed for 
suppressing the oscillating electrostatic field without de- 
stroying the DL structure. The proposed mechanism is 
well confirmed by Particle-in-Cell (PIC) simulations, and 
the carbon ion beams driven by the elliptically polarized 
lasers own much better energy spectrum than those by 
circularly polarized lasers. 

Fig. 1 (a) illustrates the classical HBRPA schem o 19 ' 20 . 
The ponderomotive force of a laser pulse drives the local 
electrons inward, resulting in a shock-like DL region with 
huge electrostatic charge-separation field. This electro- 
static field traps and reflects the ions ahead of the DL, 
and compresses and accelerates them inward like a pis- 
ton. We consider in the piston-rest fram o 19 ' 20 . The ions 
which incident toward the DL with the velocity (3f are 
reflected at x' = —D[, where /?/ is the piston moving 
velocity and D[ is the width of the DL in the piston-rest 
frame. The electrons and laser light arc reflected at the 
x' = 0. From the momentum balance we have 19 
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FIG. 1. (color online) Schematic structure of the electrostatic 
DL maintained by the radiation pressure in the piston-rest 
frame, for (a) CP laser and (b) EP laser. 




FIG. 2. (color online) The maximum electrostatic DL field 
Ex {%' = 0) vs. time, where the thin black (red) curve is for 
CP (EP) laser pulses. The thick black (red) line shows the 
maximum electrostatic DL field in the absence of the high- 
frequency components of the CP (EP) lasers. The laser in- 
tensity is 4.4 x 10 21 W/cm 2 with 1.0 /im wavelength, carbon 
plasma density is n e = 20n c and the ion charge number is 
Z = 6. For EP laser pulse, the polarization ratio is a = 0.64. 



where mj, m e , rii, and Z are the ion mass, electron mass, 
ion density, and ion charge number, respectively, 7/ = 

s 3 is the CP laser intensity, and 



In = aim e n r 



ao = eEo/m e u}QC is the normalized CP laser amplitude. 
The velocity of the piston normalized by the light speed 
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FIG. 3. (color online) (a) The maximum electrostatic DL field 
and width according to Eq. (6) and (7) vs. laser intensity, here 
the DL width has already transformed to the one in laboratory 
frame, (b) The optimum polarization ratio a = a z /a y vs. the 
laser intensity. The carbon plasma density n e = 20n c and ion 
charge number Z — 6 is given. 



c can then be written as j5f — /lnjrii(mi + Zm e )c 3 / 
[1 + / Iq / rii{mi + Zm e )c 3 ] . As the ion acceleration takes 
place mainly in the DL region, the equation of motion of 
an ion particle within the DL is 20 

m^d/Jdx' = ZeE x , (2) 



where 7- = — (i'^ and E x is the electrostatic DL 

field. The Poisson equation is 

dE x /dx' = AitZen[. (3) 

Ion continuity leads to v! i — 2riij f f3 f / f3[. From Eqs. (2) 
and (3), we can then obtain 

d/'Jdx' 2 = 2 Wp 2 i7/ /Vc 2 #. 



(4) 



At the reflection point x' = —DL we have fi[ = 0, so that 
7- = 1. Charge conservation leads to E x (x = —D[) = 
or d/Jdx' = 0. 

The first integral of Eq. (4) after multiplying dj/dx' 
on both sides is 



d/Jdx' = 2u pl /Wfil? ~ l) V4 /c, 



(5) 
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which relates the ion kinetic energy to the electrostatic 
DL field. At the position x' — of peak electron density, 
7- can be approximated by 7-(x' = 0) = 7/ since E x 
outside the DL region is very small and has little effect 
on the incident velocity /3f. Accordingly, the maximum 
electrostatic DL field E x (at x' — 0) is 

E xmax = {m, l c 2 /Ze)d'j'Jdx' = 2rrncuj pl f3f"ff /Ze. (6) 

Integrating Eq. (5) and assuming /3f <C 1, we can obtain 
the double layer width approximately as 20 

D{ = I3 f c/3uj pi j f , (7) 

which corresponds to Di — "ffD' t — /3fc/3uj p i in the lab- 
oratory frame. 

The analysis above assumes quasistationarity in the 
piston-rest frame. In reality, as can be seen in Fig. 2, the 
maximum electrostatic DL field E x is oscillating 20,21 - 25 . 
In fact, in the piston-rest frame E x does not instanta- 
neously reflect the ions. Instead, ion acceleration involves 
a finite tim o 20 ' 21 . When the ions stream into the DL 
region, the total charge density grows and a large elec- 
trostatic field is induced. However, when the ions are 
reflected, the total charge density and the electrostatic 
field decrease. Thus the sawtooth structure of the elec- 
trostatic field is closely associated with the steep varia- 
tion of the ion density in the DL region. The large os- 
cillating electrostatic DL field disperses the accelerated 
ions 2 ^, which would otherwise be mono-energetic. 

One can expect that if the reflections of ions and elec- 
trons are synchronized, as depicted in Fig. 1 (b), the 
sharp variations of the ion density and oscillations of the 
electrostatic DL field can be reduced or eliminated. Here 



m( 7/ - 1)/Z + yjl + (a 2 - a 2 ) 2 /16 - 1 = 



we propose to accomplish this by using an EP laser pulse. 
In this case, the laser intensity in Eq. (1) is replaced by 
I Q = (a 2 + a 2 z )m e n c c 3 /2. The J x B effect of the EP 
laser pulse drives the electrons out of the x' = plane, 
and fills the double layer region with hot electrons. The 
longitudinal electron momentum is 2 ^ p e — (a 2 — a 2 .)/ 4, 
with a y > a z . The piston structure can be maintained 
if these electrons are stopped (in the forward direction) 
within the DL. Balancing the (forward) kinetic and the 
electrostatic and ponderomotive potential energies of the 
electrons, we can obtain the most suitable polarization 
ratio a. The ponderomotive force acting on the electrons 
is^ 

2(Vc)(l - /?/)/(! + Pf) = J f P <(x')dx', (8) 

where the integral can be approximated by f p n e c/oj pe . 
Thus, the averaged ponderomotive force f p is 

f p = 2u pe {I Q /n e c 2 ){\ - p f )/{\ + p f ). (9) 

From the energy balance relation, we can get the opti- 
mum coupling condition 



J eE x (x')dx' + m e c 2 ( le -l) = f p D' h (10) 
- D l 

where 7 e = y/l + p 2 . According to Eq. (2), the integral 
can be approximated by rriiC 2 { r )f — 1)/Z if the effect of 
the very weak E x outside the DL region is neglected. The 
condition Eq. (10) can then be written as 



V^/Z(al + a 2 z )(3 f (l - (3 f )/3(l + p f ) 7f n, (11) 



where m = mi/m e and n — n e /n c . For given plasma den- 
sity and laser intensity which is Iq — (a 2 + a 2 )m e n c c 3 /2, 
the optimum polarization ratio a (which equals a z /a y ) 
for the proposed scheme can then be obtained from Eq. 
(11). Fig. 3 (b) shows the optimum polarization ratio 
a as a function of the laser intensity for given carbon 
plasmas with density n e — 20n c and ion charge number 
Z = 6. 

We have applied a series of one dimensional (ID) 
PIC (KLAP cod o 30 i 31 ) simulations to further confirm this 
mechanism. As for the longitudinal problem, ID PIC 
is enough to demonstrate the main physics. The size of 
the simulation is 8OA0 with Ao represents the laser wave- 
length. The simulation box is divided into uniform grid 
of 16000. The CP laser pulse propagates into the sim- 
ulation box from the left boundary. The thick target 
consists of two species: electrons and carbon ions with 
ion charge number z = 6, which are initially located in 



the region 10Ao < x < 40Ao with density n e = 20n c , 
where n c = WQe 2 m e /47r = 1.1 x 10 21 /cm 3 is the critical 
density for 1.0 /im wavelength laser pulses. We use 600 
electrons and 100 carbon ions per cell to run the simula- 
tions. The normalized amplitude of the CP laser electric 
field is a y — a z — 40, corresponding to the laser intensity 
4.4 x 10 21 W/cm 2 . The CP laser pulse is of constant tem- 
poral intensity profile. In contrast, the EP laser pulses 
with the same intensity and temporal profile are also run. 
Here in our simulations, we have scan the EP polarization 
ratio a to obtain the most suitable one with the narrow- 
est energy spectrum, which is a = 0.64. While according 
to Eq. (11), for the given laser intensity and carbon ion 
density, the theoretical predicted one is a = 0.68, which 
is quite close to the PIC optimized one. 

Fig. 4 (a) and (b) show the profiles of the electron, 
carbon ion densities and the electrostatic DL field driven 
by the CP and EP laser pulses, respectively. We see that 
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FIG. 4. (color online) The electron, carbon ion density and 
electrostatic DL field profiles corresponding to the CP and 
EP laser pulses, respectively, at t — 50To . The simulation pa- 
rameters are the same as for Fig. 2. The insets are magnified 
views of the DL regions. 



these profiles are consistent with the physical pictures as 
shown in Fig. 1. The effective ion-acceleration electro- 
static field corresponding to the CP laser pulse is mostly 
within the DL region. From our simulation the maxi- 
mum electrostatic DL field and width for CP laser pulse 
are about 60.0 normalized unit and O.lOAo which agree 
well with the theoretical predicted ones as shown in Fig. 
3 (a). As discussed, the electrostatic DL field is sawtooth 
oscillating. On the other hand, when the EP laser pulse 
is used, the reflection of electrons and ions take place 
together. In this case the width of the DL becomes prac- 
tically zero. However, the space-charge electric field is 
nonzero since the electron density at the reflection point 
(region) remains lower than the ion density. In fact, it 
has a long tail in the laser propagation direction as the 
ions and electrons move together. We can see in Fig. 2 
that the ion-accelerating electrostatic DL field is indeed 
not oscillating, which is the key to realizing the mono- 
energetic ion beams. 

Here we go to the details of the oscillating electrostatic 
DL field. Fig. 2 (c) shows the time evolution of the max- 
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FIG. 5. (color online) Energy spectrum the carbon ion beams 
from PIC simulations, (a) The energy spectrum of proton 
beams driven by CP and EP laser pulses at t — 40To , t = 60Tb 
and t = 80Tb, in which the polarization ratio is a = 0.64. (b) 
The energy spectrum of the carbon ion beams driven by EP 
laser pulses with polarization ratio a — 0.61, 0.63, 0.64, 0.65, 
0.67 at t = 80T . Here, the laser intensity is 4.4 x 10 21 W/cm 2 
with 1.0 /im wavelength, carbon plasma density is n e = 20n c 
and the ion charge number is Z = 6. 



imum DL electric field for both the CP (black curve) 
and EP (red curve) laser pulses. Two oscillation compo- 
nents can be identified. The high-frequency component 
is induced by the electron motion, and the low-frequency 
oscillations at several laser periods are induced by the 
trapping and reflection of the ions^i Because of their 
large mass, the carbon ions cannot affect nor respond to 
the high-frequency component of the electric field. That 
is, the ions respond only to the time averaged (over the 
fast variations) electric field, as shown by the thick black 
and red curves. We can see that the sawtooth electric 
field profile corresponding to the CP laser remains, but 
that corresponding to the EP laser is greatly suppressed. 

As expected, after suppressing the longitudinal elec- 
trostatic DL field using EP laser pulses, the carbon ion 
beam quality is greatly improved compared with that of 
CP laser pulses, which is clearly shown in Fig. 5 (a). 
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For the ion beams driven by EP laser pulses, the peak 
energy can be as high as 360 MeV with half-width-half- 
maximum energy spread less than 5.0%. To test the ro- 
bust of this scheme, EP laser pulses with polarization 
ratio 0.67 > a > 0.61 have been run, the quality of the 
carbon ion beams remain much better than that driven 
by the CP laser pulse with the same intensity and tem- 
poral profile, which is shown in Fig 5 (b). 

To achieve the fast ignition (FI) driven by ion beams, 
the specification of the beam quality is as follows 3 ^— : 
10 16 proton particles with energy 3^5 MeV or car- 
bon ion 10 14 with energy 300 - 500 MeV. The relatively 
smaller particle number makes the FI driven by carbon 
ions more attractable, but it calls for higher standard of 
the energy spread. For example, if the peak energy of 
carbon ion beams is 400 MeV and the distance (from ion 
source to the compressed core) is 2 cm, 10% and 50% 
energy spreads correspond to 12.5 kJ and great than 20 
kJ ignition energy 3 ^. We believe that the EP laser driven 
HBRPA is a potential way to realize the specification of 
the beam quality for FI, which is of large particle number 
and narrow energy spread. 

It should be mentioned that the scheme has also been 
applied to proton-acceleration, and the EP laser pulses 
are confirmed to be quite efficient than the CP laser 
pulses. If the laser wavelength is longer, say 10.6 /xm 
of CO2 lasers, the laser intensity can be reduced to be- 
low 10 20 W/cm 2 , but the plasma density should be re- 
duced accordingl y 26 ' 35 . From the relations (1) and (11), 
the quality of the ion beams would not change, but the 
particle number would be significantly less. 

For thin target, usually the LS model of RPA domi- 
nates. However, the acceleration mechanisms of LSRPA 
and the HBRPA are quite different. For LSRPA, the elec- 
trons and ions form a neutral layer, which, as a whole, 
is continually accelerated by the light pressure. Once 
the neutrality of the layer is ruined, the ion particles 
would undergo Coulomb explorations, and the acceler- 
ation would terminate. The CP laser, with constant 
ponderomotive force, can well keep the neutrality of the 
layer, which means that the CP laser is the best choice 
for LSRPA 3 ^. The energy and intensity of accelerated ion 
beams diminish with decreasing a and almost disappear 
for a < 0.7. 

In summary, we have proposed to use an EP laser pulse 
for ion acceleration in the HBRPA regime. An analyt- 
ical model is used to estimate the most suitable laser 
polarization ratio. In this scheme, J x B electron heat- 
ing causes the DL region to be filled with an appropriate 
distribution of hot electrons such that large ion density 
variations and electrostatic field oscillations in the DL re- 
gion are suppressed, so that a mono-energetic ion beam 
can be generated. This scheme is further confirmed by 
ID PIC simulations, and the ion beams driven by the EP 
lasers own much better energy spectrum than those by 
CP lasers. 
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